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Abstract 

We consider non-colliding Brownian bridges starting from two points and re- 
turning to the same position. These positions are chosen such that, in the limit of 
large number of bridges, the two families of bridges just touch each other forming a 
tacnode. We obtain the limiting process at the tacnode, the (asymmetric) tacnode 
process. It is a determinantal point process with correlation kernel given by two 
parameters: (1) the curvature's ratio A > of the limit shapes of the two families 
of bridges, (2) a parameter a G K controlling the interaction on the fluctuation 
scale. This generalizes the result for the symmetric tacnode process (A = 1 case). 

1 Introduction and results 

Systems of non-colliding Brownian motions have been much studied recently. They arise 
in random matrix theory (see e.g. [17,18,20]), as limit processes of random walk, discrete 
growth models, and random tiling problems, see e.g. [9-13, 16,21,23,24]. 

Considering non-colliding Brownian bridges (as well as discrete analogues), various 
kinds of determinantal processes appear naturally. Assume that the starting and ending 
points are chosen such that in the limit of large number of bridges occupy a region 
bordered by a deterministic limit shape (see Figure 1 for an illustration). Then, inside 
the limit shape (in the bulk) one observes the process with the sine kernel, see e.g. [21]. 
At the edge of the limit shape, the last bridge is described asymptotically by the Airy 2 
process [10,12,23]. Whenever there is a cusp in the limit shape, then the process around 
the cusp is the Pearcey process [4,7,22,27]. All these process are quite robust, in the 
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sense that by moving the initial and/or ending points of the bridges, the only changes 
are geometric (e.g., the position and direction of the edge/cusp changes and numerical 
coefficients in the scaling) but the processes are the same without free parameter. 

The case of the tacnode is more delicate and the limit process is described by two 
parameter. Recently, three different approaches have been used to unravel the tacnode 
process. In the first work, Adler, Ferrari and van Moerbeke [I] derived the symmetric 
tacnode process from a limit of non-intersecting random walks. Meanwhile two other 
groups were after a solution for the Brownian bridge setting. Soon after [I], a solu- 
tion appeared in term of a 4 x 4 Riemann-Hilbert problem by Delvaux, Kuijlaars and 
Zhang [8]. Their solution is for the generic tacnode process. The third approach, leading 
to, in our opinion, the simplest of the three formulations was posted more recently by 
Johansson [15]. In the latter the asymptotic analysis was restricted to the symmetric 
tacnode. In the present paper, we analyse the general case starting with the result on 
two sets of Brownian bridges of [15]. 

The equivalence between the last two formulations follows from the fact that the 
starting model is the same. However it seems hard to compare the analytic formulas 
directly. The equivalence of the results between the random walk and Brownian bridge 
case is expected by universality, and it can be indirectly checked by analysing a discrete 
model with the two approaches, as it was made very recently for the double Aztec 
diamond in [2]. 

Now we introduce the model and state the result of this paper. We consider (1 + X)n 
non-colliding standard Brownian motions with two starting points and two endpoints 
where A > is a fixed parameter. More precisely, n of the Brownian motions start at a\ 
at time and arrive at a\ at time 1, the remaining 1 An Brownian particles have starting 
and ending points at ci2 at time and 1 respectively with a\ < a<i- For finite times 
ti, . . . ,tk, the positions of the particles at these times form an extended determinantal 
point process (for more informations on determinantal point processes, see [3, 14, 19, 
25,26]). For a fixed integer n and a fixed A > 0, let us denote by £ n ^\ n {s,u,t,v) the 
kernel of this determinantal point process with s, t G (0,1) and u, v G R. The kernel 
t,v) was obtained in [15], see Theorem 2.1 below for the formula. 

In this paper, we take the n — > oo limit in the model described above. The global 
picture is that the two systems of non-colliding Brownian motions form two ellipses 
touching each other at a tacnode (see Figure 1 for an illustration). Under proper rescal- 
ing, we obtain a limiting determinantal point process in the neighborhood of the point 
of tangency. 

Here we consider the general case when two parameters modulate the limit process. 
One of them is the strength of interaction, called a, the other one is a measure of 
asymmetry, called A which we have chosen to be the ratio of curvatures of the two 
ellipses at the point of tangency. For A = 1, we get back to the symmetric case treated 

! We do not write here integer part of An to keep the notation simple. 
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Figure 1: The asymmetric system of non- colliding Brownian motions with 15 respec- 
tively 30 paths in the two groups, i.e. n — 15, A = 2. 



in [15]. 

The scaling of the starting and ending points is 



a 1 = -(Vn+|n- 1 /«), (1) 
a 2 = yft (yE + In- 1 '*). (2) 



We denote by a the distance of the two endpoints: 



VA) (v^+£n- 1/6 V (3) 



a = a2 — cti = 

In this setting, the tacnode is at (1/2,0), so that the space-time scaling we need to 
consider is 

s = l(l + T 1 n- 1 ' 3 ), t = i(l + r 2 n- 1 / 3 ), 

1 1 (4) 

The limiting kernel takes the form 

^tac (Tl, 6, 72, 6) = < T-2)p(r 2 - Ti; 6) 
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where 

^;x, 2 /) = - 7 L=exp(-^ 1 ^) (6) 

is the Gaussian kernel. To describe L tac , we need to introduce some notations. For a 
parameter s, let 

Ai (s) (x) = e 2 z s3+xs Ai(s 2 + x) (7) 

be the extended Airy function where Ai^ = Ai is the standard Airy function. The 
extended Airy kernel is given by 

POO 

K { ^\x,y) = / Ai {a) (x + u)Ai w (y + u)du (8) 
Jo 

where K^ ^ 1 = is the standard Airy kernel. Let us denote the function 



B^(x) = Ai (r) [i + (l + Vx^^j V3 fij Ai(x + n)dn (9) 

which is reminiscent of the definition of the Airy kernel. Let also 

b^(x) = X^A^i-X 1 ^ + (1 + v/A) 1/3 x), (10) 

and define 

i&(ti,Zi, t 2 , 6) = 47 T1,T2) (o- + 6,0- + 6) 

+ (1 + v^) 1 / 3 - fc*^, (1 - x~«K^y l B* T1>aHl ) L2{{ z i00)) (11) 

where Xa{%) — l(x > a) and 

3 : = A 1 / 6 (l + v / A) 2/3 a. (12) 
Now we can state our main result. 

Theorem 1.1. The (asymmetric) tacnode process T CT ' A obtained by the limit of the two 
non- colliding families of n respectively Xn Brownian motions under the scaling (l)-(4) 
in the neighborhood of the tacnode is given by the following gap probabilities. For any k 
and ti, . . . , tj, G (0, 1) and for any compact set E C {ti, . . . , x M., 

P(T^(1 E ) = 0) = det(l - C^) L2(E) (13) 
where C^ c is the extended kernel given by (5). 



4 



Remark 1. The (asymmetric) tacnode process has an intrinsic symmetry under the 
reflection on the horizontal axis that is inherited from the finite system of Brownian 
motions. This corresponds to the following transformation of the variables: 

A — ► A" 1 , (14) 
n — > An, (15) 

n — ► \ l ^r h (16) 



e . _^ _ A l/6 6> (17) 

a — > A 2 /V. (18) 



The different powers of A in the change of parameters (16)-(18) is necessary for observing 
the process on the same scale. Note that a given in (12) is left invariant under the above 
transformation. 

Next we present an alternative formulation of C^, inspired from the analogue re- 
formulation of the kernel in [2]. Let us introduce the function 

C* € (s) = b x ^(x) - B^(x) 

= A 1/6 Ai (AV3r )(-A 1/6 £ + (1 + V\Y/ Z x) (1Q) 



Ai (r) (£ + (1 + \^) 1/3 fx) Ai(x + /i)du 



Using this definition, we can give another expression for the kernel £^ac which is formally 
similar to (5), but the ingredients can be given by a single integral as follows. 

Proposition 1.2. With 

POO 

L t ^(T 1 ,ei,T 2 ,e 2 ) = (l + VF T ) 1 / 3 / (l-xzKMX^C^^Wlkrtetote, (20) 

J a 

we have 

^tac( r l^l,T"2,6) = ~Mn < T 2 )p(T2-T 1 ;€ 1 ,£ 2 ) 

+ L^(n, £i, r 2 , 6) + X^I^**"* (A i/3 Ti) _ A i/e eij ^ _ x i/ %) _ (21) 
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2 Johansson's formula 

In this section, we recall Theorem 1.4 of Johansson in [15], the starting point for our 
analysis. He obtains a formula for the correlation kernel of two non-colliding families of 
Brownian particles with the following properties. The first family consists of n particles, 
and they start at position a± at time and end at a\ at time 1. The other family has 
m Brownian motions which start at position a 2 at time and end at position a 2 at 
time 1 with a 2 > ai- This system of Brownian motions conditioned on no intersection 
in the time interval (0, 1) forms an extended determinantal point process with kernel 
£n,m(s,u,t,v) given as follows. 

Let a = a 2 — d\ and d > a, parameter which can be chosen freely in Theorem 2.1. 
We use the notation 

Autv= = / dw dzi- '—^-) 22 

sz 2 uz tw 2 vw 



x exp I — -— r — a\Z + h — r + a\W 



2(1-5) 1 1-5 ' 2(1-0 ' 1 l-t / 

B\ v {*) = } d f^— fdwf dz( - z/a 2 r — (23) 
MV (2m) 2 VT^tJiR J Dai \l-zfaj z-w 

( tw 2 vw 

x exp I — + diW — h axz 

V 2 f 1 t j 1 t 

_i , . dy/a f ( w\ m ( tw 2 vw \ 

Ux) = wn L dt " I 1 " ^ ) exp [in— +a ' w -— t + axw ) ' (24) 



(2m) VI - s J Dai J Da2 \l-z/a x J \l-w/a 2 J w-z 



xex P |-^ r — J -a 1 ^+ r — -ayw 

Mfo,y) = ^ / d* / d^T^yf^-^r— (26) 
01 (2m) 2 J Dai J Da2 \l-z/aj \l-w/a 2 ) z-w 1 J 

where _D ai and _D a2 are counterclockwise oriented circles around a\ and a 2 respectively 
with small radii. 

Very similarly, let 



{2 = d f dw f dz( l ~ w l a2 

s ^' v (2m) 2 ^(l-s)(l-t)J tR W J Da2 \l-z/a 2 ) w-z 

sz 2 uz tw 2 vw 
x exp I —777- r — a 2 z + h — r + a 2 w 



2(1 -s) 1-8 2(1 -t) 1-t 

(2m) 2 VT=t J lR J Da9 \l-z/a 2 ) z-w 



tw 2 vw 

' 2(1 - 1) + fl2W; ~ r^t ~ axz 

^ (x) = 2Wmi K d "l 1 "^J ex H^) +a2 "~T^- aa;w J' (29) 

CUV)= , d ^— f dzf d W (— VVf^T^r^- (3°) 

(SZ^ uz 
-— r - a 2 z + h ayw 
2(1 — s) 1 — s 

M 2 ( X ,y) = T ^-[ dzf dw(^^Y(\^pY^—e-^ ^. (31) 
Furthermore, let 

1 (in — t>) 2 u 2 v 2 \ 

i,) = vm^T) exp vw^> + - w^T)) 1(5 < () (32) 

denote the conjugated Brownian kernel. 

The theorem below is a consequence of Theorem 1.4 in [15]. 

Theorem 2.1 (Johansson 2011). The extended determinantal kernel for n + m non- 
colliding Brownian motions with two starting points and two endpoints described above 
can be written as 

£„, m (s, u, t,v) = - q{s, u, t, v) 

+ d- 2 Al >Ujt>v + d~ 2 ((Bl + (1 - MlY l Cl u ) L , {{hoo)) (33) 
+ d~ 2 A 2 M + d- 2 ((B 2 t>v + (1 - M 2 )- 1 ^)^^. 

Remark 2. It is proved in Lemma 1.2 of [15] that det(l — M^^m^ > and 
det(l — Mq) L 2(( 1oo )) > so that (33) makes sense. 

Remark 3. This configuration has a natural symmetry. By reflecting the vertical 
direction, one observes the same process with parameters modified as follows: 

n -H- m, a± — > —a 2 , a 2 —> —a\, u —> —u, v — > —v, (34) 

s and t are unchanged. It is easy to check that the ingredients A], utv ,B] v , f3} v) C] u 
and Mq of the kernel of the finite system in (22)-(26) transform to their counterparts 
with upper index 2 after taking the change of parameters (34) and that q(s,u,t,v) is 
invariant under this action. 

The symmetry of the limiting tacnode process established in Remark 1 is a conse- 
quence of the discrete symmetry. 
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3 Proof of the main results 



First, we give the proof of Theorem 1.1 using two lemmas which are proved in Section 4. 
We start with formula (33), and we apply it to our present setting. Then, we perform 
asymptotic analysis for the functions obtained in this way. 
The appropriate order of the parameter d is 

7,-1/12 

rf=- ? =-, (35) 

since we want d 2 C n> \ n to converge, so d 2 is the scaling of the space variables, see (4). 

The strategy of the proof is that first, we establish pointwise convergence of the 
elements of the kernel to the appropriate functions in Lemma 3.1. Then, we give uniform 
bounds on the functions in Lemma 3.2. This gives, using dominated convergence, that 
the kernel £ n! \ n under the scaling (l)-(4) converges pointwise to ^ turns out that 
this convergence is uniform on compact sets. Dominated convergence ensures that also 
the gap probabilities expressed by Fredholm determinants converge. 

Lemma 3.1. Under the scaling given by (l)-(4) ; the following pointwise limits hold as 
n — » oo. 

n'^Blil + xn" 2 ' 3 ) -+ -(1 + V\) 1/2 B^ H2 (a + x), (37) 

n- 1/3 /^(l + xn 2 ' 3 ) -+ (1 + V\) l,2 b x T2 ^ + (38) 
n~ 1/3 C s yi + yn^ 3 ) -+ -(1 + Vxf/ 2 B\^{a + y), (39) 



n 



~ 2 / 3 M^(l + xn- 2/3 , 1 + yn- 2 '*) ->• A 1/6 (l + VX) 2/3 K Ai {a + x,a + y). (40) 



where x = A 1/6 (l + \fX) 2/3 x and y = A 1/6 (l + \/X) 2/3 y and see (12). 
The analogue for the second set of terms is 

AL, t , v -> X^K^^iX^a - X^%, A 2 /V - A 1/6 £ 2 ), 

(41) 

^Bl(l + xn- 2 ' 3 ) -> -A 1/6 (l + VX^B^^^a + x), (42) 

+ XU- 2 / 3 ) -> A 1/6 (l + V^Xvl^/a^/e^ + 5), (43) 

- 1/3 C(! + ^~ 2/3 ) -+ - AV6 (! + ^W^av. 6 P + y), (44) 

n-^Mo^l + xn" 2 / 3 , 1 + yrT 2 / 3 ) -> A 1/6 (l + y/X) 2 / 3 K M {a + x,a + y). (45) 
TTie convergence is uniform for £i and £2 zn compact subset ofM.. 



n 
n 
n 
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\A i 1 < C 




\n- 1/3 B{ v 


;i- 


-xn~ 2/3 )\ < Ce 


- cx 




;i- 


-xn~ 2/3 )\ < Ce 


- cx 


\n~ 1/3 Clu 


;i- 


Vyn~ 2,3 )\ < Ce 


-cy 


2 / 3 M*(l + xn" 2 / 3 , 


i - 


Vyn~ 2l3 )\ < Ce 


-c(x+y) 



Lemma 3.2. There are constants c, C > such that for all x,y > 0, we have the 
following bounds 

(46) 
(47) 
(48) 
(49) 
(50) 

for i = 1,2. These bounds are uniform for £i and £2 in a compact subset ofM.. 

Proof of Theorem 1.1. First, we show that, with the scaling (l)-(4) and (35), we have 

d 2 C nM (s,u,t,v) — > £ t A a ^(T-i,£i,T2,£ 2 )- (51) 

The convergence of d 2 q(s, t, u, v) to the first term on the right-hand side of (5) is obvious. 
By (36) and (41), it is enough to work with the scalar products in (33). 
As in the original formulation in [15], we write 

det (l - Ml + (B\ v + &y ® C s 1 ) J i2((1)0o)) 



((^ + A!J,(1-K)~ 1 Ol 2 ((1,oo)) 

For this proof, let 



det(l - MoVtfi, 



00)) 



V(x,y) = M'(x,y) - (j^„(x) + ^ u (x))Cj ttt (y). 
Then the Fredholm determinant in the numerator of (52) can be expressed as 



(52) 



(53) 



y (zS 
^— ' m! 



m=0 



[l,oo) 



det (V (p u pj)) 



l<j,j'<m 



(Tp. 



This is equal to 



£ 

m=0 



ml 



det 



n 



-2/3 



xm 



-2/3 



-2/3 



(54) 

v\ 1+ ^ _ . a+ X]U s |j d m x 

(55) 

after the change of variables pj = 1 + Xin 2 / 3 A 1 / 6 (1 + -\/A) 2 ^ 3 - 

Using the pointwise convergence in (37)-(40) along with the bounds (47)-(50) and 
Hadamard bound on the determinant 2 , the dominated convergence theorem implies that 
the numerator of (52) converges to 



[0,oo) r 



As(1 + a/A)I 



£ 

m=0 



ml 



[0,00)^ 



det (D(xi, ^'))i<ij< m d m x 



(56) 



2 Hadamard bound: the absolute value of a determinant of a n x n matrix with entries of absolute 
value not exceeding 1 is bounded by 7W 2 . 
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with 



D(x, y) = K Ai (a + x,a + y) 

- (1 + v/FT 1 / 3 (B^ a+$2 (a + x)- b^ aH2 (a + x)) B^^a + y). (57) 

A similar argument can be used for the denominator of (52), which gives that the 
expression in (52) converges to the second term on the right-hand side of (11). 
In the same way one shows that 

{{%, + (1 - iKjr^Jx.M (58) 

converges to the scalar product appearing in the last term of (5) using the remaining 
set of assertions in Lemma 3.1 and also Lemma 3.2. Alternatively, one can refer to the 
symmetry established in Remark 3 to get the limit of (58). This verifies (51). 

It remains to argue that the process T a ' x exists as a determinantal point process. 
For this, we give a uniform bound on d 2 £ n) \ n (s, u, t, v) as £i and £2 are from a compact 
subset of E. 

For d 2 q(s,u,t,v), the assertion is clear, for A l sutv , it follows from (46). The two 
scalar products in (33) are bounded as follows. We again consider the right-hand side 
of (52). In the Fredholm expansion of the numerator (55), we get the bound 



-V[l + ~~ r r~r~ 7=-t, 1 + ~i~r~ 



a*(i + Va)§ I aI(i + Va)I' \*(i + V\y 



< Ce~< Xi+x ^ (59) 



with uniform C, c > on the compact subsets based on (47)-(50). Hence the convergence 
of (55) to (56) is uniform as £1 and £2 are in a compact set. The denominator of (52) is 
strictly positive by Lemma 1.2 of [15] and it converges to 

det(l - M 1 ) L 2 ((liOo)) — ► det(l - x^K AiX ^ L ^,oo)) = F 2 (a) > (60) 

where F 2 is the Tracy- Widom distribution function [28]. This shows that the first scalar 
product in (33) remains uniformly bounded on compact sets. One can proceed similarly 
with the second one. Therefore, the existence of the gap probabilities of the process T a ' x 
follows by expanding the Fredholm determinant of the finite size kernel and from the 
Hadamard bound on the determinant. This completes the proof of Theorem 1.1. □ 

The following proof is similar to that of Theorem 1.3 in [2], but the idea is adapted 
to the asymmetric case, so we give it completely. 

Proof of Proposition 1.2. In order to rewrite the kernel in (5), we define the function 



s X,c 



-(*) = AiM (A-^A 1 ^ - A 2/3 <r) + (1 + v/P^ 3 *) = X^y 3TX2/3a _ xl/ ^x) 

(61) 
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and the operator T on L 2 ((a, oo)) with kernel function 

T(x,y) = Ai(x + y-a). (62) 
One observes that since we have 

/•oo 

Kai(x, y) = Ai(cc + u — a) Ai(y + u — a) du (63) 

J a 

on L 2 ((a, oo)), one can write 

XaKMXa = T 2 (64) 

and also 

OO 

(l-X?#AiX?) -1 = X) r2r - (65) 

r=0 

Note also that, using the notations (61) and (62), we have 

cW*) = xl/6s SCl^ - TS rl(*)- ( 66 ) 

Similarly, 

One can also see easily that 

47 T1,T2 V + 6, ^ + 6) = (i + v^) 1 / 3 (s^, ^ 2 ) i2((;?)0o)) > (68) 

and also 

= (i + VX)'/» (^Cw^i^)^,, • (») 

Starting from (5) and (11), we can rewrite the kernel £ t £. using (19) as follows. 

£tec( r i,£i, r 2,6) 

= -1(ti< t 2 )p(t 2 - n; 6) 

+ Al7 n,r2 V + 6, * + 6) + a 1/6 a^7 a1/3ti ' a1/3t2) (a 2/3 o- - a 1/6 ^i, A 2 /V - a 1/6 £ 2 ) 



+ (1 + x/A^) 1 / 3 / (1 - XaKMX^C^^iC^^x) - b x _ T1 ^(x)) dx 

poo 

+ A 1/6 (l + V\) 1/3 (1 - X^AiXff)- 1 ^^^/.^/^^) 

x (C-AVa^^VScr-AVe^^) _ fr-A^n.A^a-AVe^ ( x )) 



(70) 
11 



If we use (68) and (69) for the first two Airy kernels, (66), (67) and (61) for the two 
integrals and (65), then we get 



A A ac (n, 6, t 2 , 6) = - i(n < r 2 )p(r 2 - n; f 2 ) 



*r=0 I L 2 ((a,oo)) 



, r=0 



L2((5,oo)) 



iA _1 ,A 2//3 (t 



>r=0 I L2((5,oo)) 



oo 



-(i + VA)V3/^r-«5r;;Cx./. £l ,s^ 

\r=0 / L 2 ((5,oo)) 

(71) 

Note that the scalar products in the third and the fourth terms on the right-hand side 
of (71) can be combined using (66), respectively, the second and the fifth terms can be 
joined by (67) yielding 

A A ac Oi, £i, r 2 , 6) = - l(n < 72)p(r 2 - n; 6, 6) 

+ (1 + x/A) 1 / 3 ((1 - x^aix?)- 1 ^^,^: 1 ' 



/ 3 ^,-A 1 ^ 2 / i2(( ~ oo)) 

+ (1 + v/A) 1 / 3 ((1 - X^AiX.)- 1 ^^,^^^,,, <l) L2((?oo)) 

"(72) 

which proves the proposition by (61). □ 



4 Asymptotic analysis 

Proof of Lemma 3.1. We write down the ingredients of the kernel in Theorem 2.1 with 
the values given by (l)-(2). By Remark 3, it is enough to prove the first set of state- 
ments (36)-(40). Then, we use the method of saddle point analysis to get the limits in 
Lemma 3.1. For the asymptotic analysis we use a standard pattern, explained e.g. in 
Section 6 of [5]. 



12 



For (22), after change of variables w — > aw/ (1 + vT) and z — > —azj (1 + a/A), we get 



d 2 a 



s,u,t,v 



(l + VX)y/(l-3)(l-t) (2 



7TZ 



x exp 
x exp 



a 2 sz 2 



dw / dz 

a 2 z 



1 + w 



1 



1 — z ) w + z 

auz 



(1 + \/A) 2 2(l - s) (l + v^) 2 (l + VA)(l-s 



(73) 



a 2 tw 2 



a 2 w 



(i + V\) 2 2(i-t) (i + V\) 2 (i + V\)(i-t) 

First we observe that Re(z + w) > 0, hence we can write 



z + w 



n 



1/3 



d/i. 



(74) 



By substituting this and (3)-(4) and by Taylor expansion, we obtain 

poo 

Ku,t,v =(n 2/3 + o(l)) / d/i 

Jo 



/ dwexp (n 

2?™ J iR V 



w 



x / dzexp n 

27TZ ' 



log(l + w) + — IW 

z 2 

-log(l -*) - — - z 



n 



2/ W-n 1 / 3 [a + 6 + ^ 



n 2/3 Tiz 2 - n 1/3 [a + £i + [i]z 



x exp (O (n 1/3 w 2 + w + n 1/3 z 2 + z)) . 



A change of variables w — > aw/(l + y/X) and z — > —az/(l + \/A) in (23) gives 



(75) 



do 3 ' 2 



(1 + V\)VT^t(2ri) 



x exp 



a 2 tw 2 



dw / dz 

'Di 

a 2 w 



1 + w 
1 - z 



1 + ^= 



Z \ 1 



VxJ z + » (76) 



a k 



(1 + VA) 2 2(1 - 1) (1 + v 7 !) 2 (i + Vx)(i-t) i + Vx 



which yields 



n 



POO 

1/3 Bl(l + xn~ 2 ' 3 ) = (n 2 ' 3 + o(l))(l + V\) l/2 / d/i 



dwexp 



log(l + iw) + — — w 



n 2 ^ 3 T2W 2 — n' 



1/! V + 6 + 

x — / dzexp (n A log (1 + -A=] - log(l — 2) — (1 + y/\)z ) 
x exp (-n 1/3 [(l + >/A)(sc + cr) + fi]z + O (n 1/3 w 2 + w + n 1/3 z 2 + z)j . 



(77) 
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In the definition of /3} v , after w — > —aw /{I + a/A), one obtains 



da 3 / 2 



(1 + V / A)v / T^ 27ri 
x exp 



dw 1 



w \ 



An 



a 2 tw 2 



+ 



a 2 w 



+ 



avw 



a 2 xw \ 



.(1 + VA) 2 2(l-t) (l + v 7 !) 2 (l + VA)(l-t) l + v^y 

Hence, 

n~ 1 / 3 /3 t 1 , ) (l + xn~ 

2/3 ) = (n 2/3 + o(1))(1 + v ^ ) l/2 



I • 

(78) 



x 



2tu 



du>exp I n 



A log 1 + 



w \ w 



w 



(79) 



x exp I n 2/3 T 2 w 2 - n 1/3 



[VXa - 6 + (1 + VX)x]w + O (n 1/3 w 2 + w)^j 



Similarly in (25) with w — > aw/(l + V~X) and z —> —azj (1 + a/A), we get 

An 



cUv) 



da 3 ' 2 



1 



(1 + V\) VT^ (2ttz) 2 J Dl 
a 2 sz 2 

x exp 



dz / dw 



a 2 z 



1 + w 

1 - z 



1 H_ 



w + z 
a 2 yw 

(1 + V\) 2 2(l - s) ~ (1 + V\) 2 ~ (l + V\)(l-s) ~ 1 + y/X 
After substituting (3)-(4), it becomes 

POO 

n-^Cl^l + yn^ 3 ) = -(n 2 ' 3 + o(l))(l + V\) l/2 \ dp 
1 



(80) 



x 



2ni 



dzexp I n 



ih 



log(l -z)- — -z 



n 2/3 Tiz 2 - n 1/3 [a + ^ + fj]z 



x 



/ dw exp in log(l + w) — A log (l = ) — (1 + \/\)w ) 

271 *Jd^ V L V vAy J/ 

1/3 [(1 + VX)(o- + y) + n]w + O (n 1/; V + w + n 1/3 ^ 2 + 2)) . 



x exp I — n 



(81) 

Using the same change of variables w — > aw /{I + \/A) and z — > —az/(l + vA) in 



Mq, we have 



Ml(x,y) 



dz I dw 

(1 + y/\) (2vrz) 2 y Dl y D 



x exp 



a 2 xz a 2 yw 



1 + w 
1 - 2 



An 



7a 1 



I j 

VxJ 



z + w 



(82) 



1 + Va 1 + vV ' 
11 



That is, 

-^M^l + xn 2 l\ 1 + yn- 2 '*) = (n 2 / 3 + o(l))(l + y/X) 
1 



n 



d/i 



/ du>exp n 

2m 1 



x exp 
1 

x 

2m 

x exp 



v'A 



-n 



log 1 



Alog(l--^L) -(1 + VX) 



1 / a [(l + VA)(^ + y) + A *]i 
y dzexp^n A log ^1 + -^=j - log(l - z) - (1 + V^z J 
(-n x / 3 [(l + VA)((7 + z) + //]* + O (n 1 / 3 ^ 2 + w + n^V + 2)) . 



(83) 



We can take the integration paths Di and -D^ to denote here circles around 1 and y/X 
with radii 1 and y/X, in this way, passing through 0. 
In the above formulae, integrals of form 

1 



2m 



dze nfo{z)+n fl{z)+n h{z) 



(84) 



appear. Hence, we can follow the steps of the proof of Lemma 6. 1 in [ ] . It can be checked 
that in all the cases in (75), (77), (79), (81) and (83), w c = z c = is a double critical 
point for the functions below that appear as /q. We also give the Taylor expansion for 
later use 



1 1 

TV TV 

log(l + w) + — -w = — + (w 4 ) 



w \ 



A log H — -= + — — w 



1 w 3 



VxJ 



Vx 3 



0(w 4 ) 



(85) 



logfl 



w) 



A log 



Vx) 



w 



y/Xw 3 



Vx 



and their versions after taking the transformation —fo(—w) which preserves the leading 
terms on the right-hand side of (85). 

The issue of finding steep descent paths is treated in a more general setting in the 
proofs of Lemma 3.2 and 4.1 below. Thus the details are omitted here. As it turns out, 
the integration paths in (75), (77), (79), (81) and (83) are steep descent for f in the 
following sense (7 means the generic integration path): 

• Re(/o(z)) on 7 reaches its maximum at 0, 

• Re(/o(^)) is monotone along 7 except at its maximum point and, if 7 = D\ or 
Dr^, then also in the antipodal points 2 or 2vA respectively (where the real part 
reaches its minimum). 
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In the neighborhood of 0, we can slightly modify the paths iR, D\ and in such 
a way that they are still steep descent and that, close to 0, the descent is the steepest. 
The modification is as follows. For a small 5 > that is given precisely later, we choose 
the path {e~ i7r/3 t, < t < 5} U {e in/3 t, < t < 5} close to 0. This is clearly locally 
steepest descent for all functions in (85). In the case of iR, we continue with 5/2 + iR 
to infinity. If we had a circle, we decrease the radius slightly, and we cut off a piece of it 
close to in such a way that it matches the steepest descent path. Computations which 
are done in the proofs of Lemma 3.2 and 4.1 show that these modified paths are still 
steep descent for f in all the cases in (75), (77), (79), (81) and (83). 

Along the integration paths, only the contribution on 7 D {\w\ < 5} is considered, 
since, by the steep descent property of the paths, the error that is made by neglecting 
the rest of the path is 0(exp(— c5 3 n)) as n — > 00 uniformly as £1 and £2 are in bounded 
intervals. 

We choose now 5 such that the error terms in (75), (77), (79), (81) and (83) are 
small enough in the 5 neighborhood of 0. That is, for the difference of the integrals 
with and without the error in the exponent, we apply \e x — 1| < |x|e' x '. We do a change 
of variables n l ^w — > w, and, by taking 5 small enough, we see that the difference is 
(^(ri" 1 / 3 ) uniformly. 

Finally, we can extend the steepest descent path to e~ l7r ' 3 oo and e 47r / 3 oo on the price 
of a uniformly 0(exp(—c5 3 n)) error again. Then, we apply the formula 

1 f W3 °° f Z " u 2 V -1/3 b C \ ^ . / b 2 C \ 

— exp [a— + bz + cz ) dz = a 1 exp — Ai — - - — - 

27ny e -i*/3 00 V 3 / \3a 2 a J \a 4 / 3 a 1 / 3 J (86) 

= a- 1 / 3 Ai( a " 2/3fe ) / 



aV3 

and we exactly get the limits (36)-(40). From this, along with the use of symmetry 
given in Remark 3, (41)-(45) follow immediately. All the error terms can be bounded 
uniformly in £ x and £ 2 if they are in a compact interval, hence Lemma 3.1 is proved. □ 

First, we give the following lemma with its full proof. In the proof of Lemma 3.2, we 
will use the assertion, and some similar statements which will not be spelled out later, 
because they can be shown as Lemma 4.1. 



Lemma 4.1. There are constants C, c > such that 

n l/3 



du>exp n 



2iri 

for s > large enough 



2 

w 



log(l + w) + — w 



Kin 2lz w 2 - sn 1/3 w 



< Ce~ cs (87) 



Proof of Lemma 4-1- We follow the lines of the proof of Proposition 5.3 in [6]. Since we 
are interested in large values of s, we take 

s = n - 2/3 s, (88) 
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and we define 



w 



fo(w) = log(l + w) + — - w - sw. 



(89) 



For small values of s, this function has two critical points at is 1 / 2 at first order, and 
we will pass through the positive one. Hence define 



a 



s 1 / 2 if?<e, 
e 1 ' 2 ifs>e, 



(90) 



for some small e > to be chosen later and consider the path T = a + iWL. By the 
Cauchy theorem, the integral in (87) does not change if we modify the integration path 

to r. 

The path T is steep descent for the function Re(/o(iu)), since 



— Re(/ (a + **)) = -* (l 



1 



(l + a) 2 + t 2 J • 



(91) 



Define 



>0 for a,t>0. 



Q(ot) = exp ( Re ( nf (a) + r?' Z K\a. 



(92) 



Let Ts = {a + it, \t\ < 5}. By the steep descent property of T, the contribution of the 
integral over r \ Ts in (87) is bounded by Q(a)0(e~ cn ) where c > does not depend on 
n. The integral on 1^ can be bounded by 



Q{a) 



n 



1/3 



2ivi 



dwexp (n(fo(w) — fo(a)) + n 2 ^ 3 k,i(w 2 — a z 



By series expansion, 



with 



Re(/ (a + it) - f (a)) = -jt 2 (l + 0(t)) 



7 



1 + a 



(93) 

(94) 
(95) 



After a change of variable w = a + it, (93) is written as 

,1/3 



7,1/3 r° 

/ dtexp (-7t 2 r2(l + 0(t)) (l + O (n~ 1/3 ))) 
^t 2 n 



1/3 



< Q(a 



2tt 



dt exp 



(96) 



<Q(a) 
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for 5 small enough and n large enough. The estimate above is the largest if 7 is small. 
Note that, by (95), (90) and (88), 

7^3 „ an V3 „ 5l/2 n l/3 „ s l/2 (97) 

which is large if s is large enough, so the integral in (93) is at most constant times Q(ct). 

Hence, it remains to bound Q(ot) exponentially in s. For this end, we use the Taylor 
expansion 

'w 3 „ \ 



If s < e, then 



fo(w)=( — -swj(l + 0(w)). (98) 



Q(a) = exp f(-f ru 8 ' 2 + Kl n 2 / 3 ^ (1 + 



(99) 

exp ^--,3/2 

where the first term in the exponent dominates as s is large, so this is even stronger 
than what we had to prove. 
If s > £, then 

Q(a) = exp ((ny/I (| - 5) + Kl n 2 / 3 e) (l + O (v^))) (100) 

where e/3 — s < — Hence, the first term in the exponent is about — | y/en 1 ^ 3 s which 
dominates the second term that is of order en 2 / 3 ~ s. Therefore, for a given e > 0, n 
can be chosen so large that 

Q(a) < exp \-^n 1/3 s\ (101) 

which finishes the proof. □ 

Proof of Lemma 3.2. By Remark 3, it is enough to prove all the bounds for i — 1. The 
assertion (46) can be shown as follows. Lemma 4.1 with s = a + £2 + yU applies for the 
integral with respect to w in (75) and provides an exponentially decaying bound in /1. 
If we prove a similar statement for the integral with respect to z, then (46) follows for 
2 = 1 along with the uniformity assertion for ^ and £2- We omit the details of the proof 
of the bound on the ^-integral here, because they are very similar to that of Lemma 4.1, 
but we give that, for the second integral in (75). We consider the function 

gf(z) = - log(l - z )-£-z-sz (102) 
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and the path {1 — pe*^} with < p < 1. It is steep descent because 

4r Re(#f(l - pe^)) = -p sin 0(2(1 -pcos0) + 5) . (103) 
d0 * „ ' 

>o 

For (47), we use Lemma 4.1 and the fact that, for the function 



A log (\ + -^=J - log(l - z ) - (1 + VX)z - sz, (104) 



the path {1 — pe 1 ^} with < p < 1 is steep descent. Indeed 



_d_ 

d0 



Re(rf(l - pP)) = -psin0 ((1 + VA) (l - ^ + *_ ) + 7) 



(105) 



where the last factor between the outhermost parenthesis is certainly positive. 
For proving (48), we take the function 



fg(w) = A log (l + ^ + ^ - w - sw. (106) 



The steep descent path is a + iR for a > by 



^fHo + U)) = -t(l- X ). (107) 



~-N^ — 



Finally, for the function 



fC(w) = log(l + w) - A log f 1 - ^ - (1 + y/X) 



w-sw, (108) 

we choose the steep descent path {vA — pe %<t> } with < p < \/\. The steep descent 
property is shown by 



±Re(f^V\-pe^)) = -psin0 ((1 + V\) (l - - 1 



pe^| 2 



(109) 



where the factor between the outhermost parenthesis is positive. This fact, together 
with the steep descent property of g$ (z) in (102), yields (49). 

Using the steep descent paths for fo(w) in (108) and for g® (z) in (104), we can mimic 
the proof of Lemma 4.1 to get exponential bounds on the integrals in (83). It completes 
the proof of Lemma 3.2. □ 
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